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Abstract 

We study the implications of minimal flavor violating low energy supersymmetry 
scenarios for the search of new physics in the B and Higgs sectors at the Tevatron 
collider and the LHC. We show that the already stringent Tevatron bound on the 
decay rate Bg — > fJ-^^~ sets strong constraints on the possibility of generating large 
corrections to the mass difference AMg of the Bg eigenstates. We also show that 
the Bg fJ-^l^^ bound together with the constraint on the branching ratio of the 
rare decay b ^ has strong implications for the search of light, non-standard Higgs 
bosons at hadron colliders. In doing this, we demonstrate that the former expressions 
derived for the analysis of the double penguin contributions in the Kaon sector need 
to be corrected by additional terms for a realistic analysis of these effects. We also 
study a specific non-minimal flavor violating scenario, where there are flavor changing 
gluino-squark-quark interactions, governed by the CKM matrix elements, and show 
that the B and Higgs physics constraints are similar to the ones in the minimal flavor 
violating case. Finally we show that, in scenarios like electroweak baryogenesis which 
have light stops and charginos, there may be enhanced effects on the B and K mixing 
parameters, without any significant effect on the rate of Bg ^jl'^ ^~ . 



1 Introduction 



The standard model (SM) provides an accurate description of all the results from high energy 
physics experiments, in particular precision electroweak measurements and flavor physics 
observables. These experiments put strong constraints on extensions of the SM that have 
tree-level flavor changing neutral current effects or large custodial symmetry breaking effects. 
For renormalizable, weakly interacting theories, where the new exotic particles acquire large 
gauge invariant masses so that they decouple from the low energy effective theory, these 
constraints can be avoided. Low energy supersymmetry PJ |2] is a particularly attractive 
example of this kind of theory. The minimal supersymmetric extension of the Standard 
Model or MSSM (with gauge invariant SUSY breaking masses of the order of 1 TeV) predicts 
an extended Higgs sector with a light SM-like Higgs boson of mass lower than 135 GeV [3]- 
[T3] that agrees well with precision electroweak measurements. 

However the structure of supersymmetry breaking parameters is not well defined. If 
there are no tree-level flavor changing transitions in any gauge or super-gauge interaction, 
then the deviations from SM predictions are naturally small. Such small deviations can be 
achieved if the quark and squark mass matrices are block diagonalizable in the same basis. 
For instance, this happens when the squark and slepton supersymmetry breaking masses are 
flavor independent. For these kinds of models, all flavor violating effects are induced at the 
loop-level and are governed by the CKM matrix elements, as in the SM. Many studies have 
concentrated on the properties of these minimal flavor violating scenarios (see, for example, 
Refs. 

In this article we shall analyze their flavor violating effects in two quite generic cases. In 
the first case, we consider a low energy effective theory in which the quark and squark mass 
matrices are aligned in flavor space and can be simultaneously diagonalized in blocks, as 
described in the next section. We will remain agnostic about how this effective low energy 
theory is UV completed. However, since the Yukawa-induced radiative corrections to the 
soft supersymmetry breaking parameters tend to destroy the alignment of the squark and 
quark mass matrices, this situation may be only naturally realized in models of low energy 
supersymmetry breaking, where these radiative corrections are small. We call this low energy 
scenario Minimal Flavor Violation. 

In order to study the possible effect of Yukawa dependent radiative corrections we study a 
second case, in which we assume a departure from the alignment condition by the presence of 
flavor violating effects proportional to the CKM matrix elements. These effects are induced 
by corrections to the left-handed down squark mass matrices proportional to the product 
of the up-quark Yukawa matrix and its hermitian conjugate (or, in general, powers of this 
product). We furthermore assume that the right-handed down squark masses are flavor 
independent. As we will discuss in more detail in the next section, these conditions at low 
energies are achieved, for instance, by Yukawa dependent radiative corrections, if one starts 
from flavor independent squark masses at a high energy scale at moderate values of tan/5. 
One characteristics of this second scenario is that there are flavor violating down-squark- 
gluino vertices at tree-level. Since all flavor violating effects are governed by the CKM 
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matrix elements, this scenario would also enter within the general definition of minimal 
fiavor violating models given in Ref . [21] . However, due to the presence of flavor violating 
couplings at tree level, we will denote it as non-minimal flavor violation in order distinguish 
it from the flrst scenario of flavor alignment at the weak scale, in which such tree-level effects 
are absent. As we will show, the phenomenological predictions in this scenario are similar 
to those of the flavor alignment case, unless the left-handed squarks and the gluino are very 
light. 

Apart from the structure of supersymmetry breaking parameters, the phases associated 
with them are also important. In minimal flavor violating schemes there are at least two 
phases that cannot be absorbed by redeflning the low energy flelds. For real values of the 
fi parameter, these phases can be associated with a universal phase for the gaugino masses 
and the trilinear mass parameter. In general, however, one can choose independent phases 
for the different gaugino masses and trilinear mass parameters. CP-violating phases beyond 
the CKM one are required, for instance, in models of electroweak baryogenesis [211^1211 • In 
this scenario, there could be signiflcant effects on AMg, BlZ{Bs /W^/W^) and ex because of 
the presence of a light stop and extra phases in the chargino, neutralino and gluino sectors. 
We shall comment on the effects of these new CP violating phases below. 

In this paper we attempt to develop a systematic method of treating the extra sources 
of flavor violation in the minimal and non-minimal flavor violating models described above. 
We show that the usual approach of calculating tan/5 enhanced FCNC (Flavor Changing 
Neutral Currents) effects in the Kaon sector does not agree with the exact results one flnds 
in the limit of flavor independent masses. Thus, we develop a perturbative approach that 
leads to agreement with the exact result in this limit. 

We shall emphasize the implications of the present bounds on BlZ{Bs 1^^ 1^^ ) for future 
measurements at the Tevatron collider, both in Higgs as well as in B-physics. In particular, 
we shall show that the present bound on BTZ{Bs /U^/W ) leads to strong constraints 
on possible corrections to both AM^ and the Kaon mixing parameters in minimal flavor 
violating schemes. Moreover, we shall show that this bound, together with the constraint 
implied by the measurement of BTZ{h — > 37) leads to limits on the possibility detecting 
light, non-standard Higgs bosons in the MSSM at the Tevatron collider. Throughout the 
paper we always take real values of /iA^, and therefore the Higgs sector is approximately 
CP- invariant and will be treated as such. 

This article is organized as follows. In section 2, we deflne our theoretical setup, giving 
the basic expressions necessary for the analysis of the flavor violating effects at large values 
of tan/3. In particular, we show how the flrst order perturbative expressions in the CKM 
matrix elements are inappropriate to deflne the corrections in the Kaon sector where higher 
order effects need to be considered. In section 3 we show the implications of the constraint 
on BTZ{Bs — > fJ^'^fi~) for the mixing parameters of the Kaon and B sectors in the large tan/? 
regime. In section 4, we explain the implications for Higgs searches at the Tevatron. We 
reserve section 5 for our conclusions and some technical details for the appendices. 



2 



4): 



0* 



^2 dR 



Figure 1: SUSY radiative corrections to the self-energies of the d-quarks in the mass insertion 
approximation 

2 Theoretical Setup 

2.1 The resummed effective Lagrangian and the sparticle spec- 
trum 

The importance of large tan /3 FCNC effects in supersymmetry has been known for sometime. 
The finite pieces of the one-loop self energy diagrams lead to an effective lagrangian for the 
quark-Higgs sector, valid at energy scales lower than the heavy squark masses, which has 
the generic form [13-120!, HaEI] 



(1) 
(2) 



- C^jf = 4Yd[$S* + K" (eo + 6VYt Y„)]4 + $^?^Y„< + h.c. 
-JZ^ass = ^4Yd[l + tan/3 (eo + eyYt Y,)]4 + ^^/?jY,n° + h.c. 

in an arbitrary basis. The and Iy matrices correspond to radiative contributions 
coming from the loops shown in Fig. Their exact dependence on the supersymmetric 
mass parameters is given in Appendix IA.2I 

The flavor structure of the loop correction factors are independent of their momentum 
integrations. Therefore, in an arbitrary basis, the flavor dependence of the loop correction 
parameters are the same as that of the mass matrices and Yukawa couplings. Thus, the loop 
correction factors have the following flavor structure 



Yd6VYtY, (X YdM-^YtM-^Y 



(3) 
(4) 



where matrices are the non-diagonal inverse squark mass squared matrices. Thus the 
sparticle spectrum is intimately connected to the e parameters which in turn affect the 
FCNC's. We look at two possible choices that connect the quark mass eigenstate basis to 
that of the squarks. 
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2.1.1 Minimal Flavor Violation 



This scenario is similar to tliat discussed in Refs. ^3 QHl ^1 ^] , where one assumes an 
alignment of the quark and squark mass matrices in flavor space. Therefore, in the low 
energy effective theory, the diagonalization of the quark mass matrices leads to squark mass 
matrices that are block diagonal. Using the following transformation matrices 



l = v2uL, 4 = U?Vocii, u% = V^un, d% = VUR (5) 



u 

to rotate the original quark supermultiplets into a basis where the tree level Yukawa couplings 
are diagonal, we get 

Yd = U^^YdU^^Vo; 
Y, = U^^Y.U^?; 
M-^ = V^M-^Hji; = Viug^M-^^UgVo; 

eo oc UgVoMr^Mr^V^U?; eo = UgVoeoV^Ug^ 

ey oc UgM-^M-^Ug^ ey = UgeyU?; (6) 

where the un-hatted mass and Yukawa matrices are diagonal and Vq is the tree level CKM 
matrix. Under this transformation the effective mass lagrangian becomes 

- -^^mass = ^dRYd[l + tan/5 (^eo + VieY|Yu|^Vo)]rfL + ^urY^ul + h.c. (7) 

where the eo and ev terms, defined in Eq. (jHl) (see also Appendix A, Eq. pH7j) and Eq. pH8|l ). 
are diagonal. Therefore the quark mass matrices receive off-diagonal terms proportional to 
ey at the 1-loop level and so need to be rediagonalized perturbatively. This procedure has 
been performed in Refs. pnil9j . However, the calculation of the (2, 1) and (1, 2) components 
of the neutral-Higgs-quark-quark coupling are affected by additional corrections not included 
in Refs. [17111^. In Appendix I A. II we calculate the corrected couplings which we present 
here. Defining the down-quark neutral Higgs interaction Lagrangian to be 

-C = di{Xl,j)''dic^s + h.c., (8) 
we find that the neutral Higgs flavor changing coupling, with I ^ J , takes the form 

^^^^^ " Vd{l + ef[tan/3)(l + eg tan/?) ^^^^ ^^^^ 

where we have ignored the small effects proportional to the first and second generation 
Yukawa couplings to find ej = + S^jeyyt, x^ and xf are the Higgs scalar components on 
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the neutral $°* and fields (see Appendix A, Eq. fjll7|) ) and 



■p3/ 
■n21 



r 



12 



ey(l + 63 tan/3) — €^(€3 — ej) tan /5 
1 + eg* tan (3 

^ , ^ 3, ^12 [(l + eotan/?)|l + e3tan/?|2- 

'1 + e2tan/5)|l + egtanpp 

■YVt tan/?(l + £3 tan + €2 tan /?) — eyy^ tan/?(l + 62 tan /?)^] 

{(1 + eo tan /3) 11 + 63 tan/3p- 

€2 tan P) — eyyl tan/?(l + 62 tan /?)(! + ei 
(e^)Vtan^/j |6y|\nan^/3 
1 + e^tan/5)(l + e^tan/?) jl + estan/^P 



(1 + 62 tan /3) 1 1 + e|]tan/3|2 
ey?/4^tan/?(l + e3tan/3)(l + 
ey tan /3 



ei - €2 



1 + 62 tan /? 



(10) 
(11) 

(12) 

tan [3) 
.(13) 



Here Vg// is the CKM matrix obtained after diagonalization of the one- loop mass matrix 
in Eq. (0). The relation between this matrix and Vq is given in the Appendix lA.ll Observe 
that in the limit of universal squark soft SUSY breaking masses the eo diagonal matrix is 
proportional to the identity and, in spite of their complicated form, all the F^'^ become 
equal to ey. The difference between the above expressions and those obtained before in the 
literature will be discussed in more detail below. 



2.1.2 Non-minimal Flavor Violation using the CKM matrix 

As explained in the introduction, we shall discuss a second scenario in which all fiavor 
violating effects are proportional to CKM matrix elements, and there are tree-level down- 
squark-gluino flavor violating vertices in the low energy effective theory. This scenario is 
similar to that discussed in Ref. [S^. For the present discussion, let us assume that we 
perform the diagonalization procedure in a single step under the transformation 

< = UgMi, 4 = UgVeffrfL, 4 = U^«i?, 4 = U^^i? (14) 

where instead of Vq the tree level CKM matrix we have Ve// the effective CKM matrix. 
This transformation leads to a diagonal quark mass matrix and a mass lagrangian of the 
form 

- C^^, = ^4ujYdU?[l + tan/3 (eo + eY|Y,|2)]V,^^rfi + ^urY^ul + h.c, (15) 

under the assumption that the matrices 

eo = U?t,^^u? 

ev = U?eVU? (16) 
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are diagonal [20] • The condition that UL^eYU^ is diagonal is the same as Eq. (jH)) in Minimal 
Flavor Violation. Thus we again need the u-squark mass matrix to be block diagonal in the u- 
quark eigenbasis. Therefore there are no flavor changing effects in the neutral up supergauge 
currents. 

However the assumption that UL^eoUL is diagonal differs from Eq. (jHl) in MFV. From 
the flavor structure of eo in Eq.Q, we see that this can only be naturally fulfilled if 

M-2 = U.^^M-^U^^, and M-^ = U^^M-^U^ (17) 

are diagonal and [Mr^, YdV^j^] = 0. The obvious way of satisfying this commutation 
relation is to require the right-handed d-squark mass matrix to be flavor independent or 
M? oc I. Observe that this analysis was not performed in Ref. [201 and hence the above 
conditions were not required in that work. As stressed in the introduction, the above flavor 
structure of mass matrices may be achieved by Yukawa induced radiative corrections to 
universal, flavor independent squark masses at high energy scales, at moderate values of 
tan/?. Assuming the squark masses are flavor independent at high energies, the only one- 
loop corrections that violate flavor are induced by the up and down Yukawa matrices because 
the gauge interactions are flavor blind. These corrections are given by [22] 



AM, 



Q 87r2 



\2ml + M^JO) + Al) Y^Y^ + {2ml + <.(0) + Al) yJy,] log f-^ 

J \MsuSY, 

(18) 

where Q denote the left-handed squarks, mo is the common squark mass at the scale M at 
which supersymmetry breaking is transmitted to the observable sector, Mfj ^(0) and Aq are 
the Higgs soft supersymmetry breaking masses and squark-Higgs trilinear mass parameters 
at that scale, and Msusy is the characteristic low energy squark mass scale. 

Similarly, the right-handed up and down squark mass matrices, receive one-loop Yukawa- 
induced corrections proportional to 

= --^ {2ml + MlS^) + Al) Y^Y^ log (j^) , (19) 

'^'^l = + + o y^yi [j^) . (20) 

respectively. 

Therefore, while the Yukawa induced radiative corrections to the right-handed squark 
mass matrices mantain the alignment of these matrices with their corresponding Yukawa 
matrices, the corrections to the left-handed squark masses induced a misalignment between 
the quark and squark mass matrices governed by CKM matrix elements. Since the dominant 
effects are governed by the third generation Yukawa eigenvalues, the down-quark Yukawa 
effects may be neglected at small or moderate values of tan/3 where the bottom Yukawa 
coupling is much smaller than the top quark one. In this case, one arrives at the properties 



and 
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of the squark mass matrices specified in the non-minimal flavor violating scenario defined 
above. 

In general, even at larger values of tan (3, the only flavor violating squark-gluino vertices 
will be in the left-handed couplings (and the Higgs-squark-squark vertices) and they will be 
governed by CKM matrix elements. The only difference between the large tan (3 case with 
respect to the non-minimal flavor violating model defined above is that the masses of the 
right-handed down squarks will no longer be flavor independent at low energies and therefore 
the eo matrix will not be aligned with the ey one. However, the flavor properties of the large 
tan f] scenario are quite similar the non-minimal flavor violating scenario specified above 
and therefore this scenario will allow us to study the possible effects of the Yukawa induced 
radiative corrections to the squark mass matrices, in particular the ones associated with the 
flavor violating down-squark-gluino couplings at tree-level. 

Following the argument in Ref. [201 '^^^ rewrite the effective lagrangian in terms of 
the mass eigenstates as 

/2 - /2 - 

-Ceff = — ($°* - tan/5)(ifimdyJ^^R" Ve//dL + — ^^^d^mdrfL 

+^IurY^ul + h.c. (21) 

where V^// is the effective CKM matrix, Yu is the diagonal up Yukawa matrix, rfid is the 
diagonal down-quark running mass matrix, and 

R= l + eotan/3 + eY|Yuptan/?. (22) 

Therefore, neglecting^ yu and Uc as compared to yt, and defining 

ej = ei + eyy^^^a ^23) 

for all J, we find 

{Yl-y = ^- -5'' (24) 

^ ^ l + ejtan/3 ^ ' 

If we assume a generational mass splitting so that the first two generations are equally 
massive and heavier than the third generation we find = eg = eo. In this case the flavor 
changing effects are not solely dependent on ey, but they also depend on the difference 
between the loop factors (ea — eo): 

(X^ = "^-^J^^g - eo)(xf - xf tan/3) gj^ gj , . 

^ .;rf(l + eotan/5)(l + e3tan/3) '^^^ ^ ^ 

The reason we call this scenario non-minimal flavor violation is that the diagonalization 
procedure induces flavor changing effects in the gluino-quark-squark couplings which lead to 

^This approximation breaks down in the limit 1 + eotan/? — > 0, the singularity in [^^f,] proportional to 
Ut cancels against those coming from j/c and y„ as discussed in Ref. |2(J| 
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additional contributions to flavor changing processes. Indeed, the assumption that eo and ev 
in Eq.()16|) are diagonal leads to the appearance of CKM elements in the down quark-squark- 
gluino coupling, as it is clear from Eqs. (fTlj) and (fT7j) . Because the left-handed squarks 
are not diagonalized by the same rotation as the left-handed quarks, the effective gluino 
Lagrangian becomes 

= -V2gs[uLrnL-URrnR] 

+ V2gs [dL r T'^ YdL - dn ^ T^dn] ■ (26) 

The appearance of the CKM matrix in the gluino couplings induces flavor changing box 
diagrams that can in principle produce large effects. 

2.1.3 The uniform squark mass limit 

The two flavor changing scenarios discussed above coincide for the case of uniform squark 
masses. Since, in this limit, the transformation performed in Section l!^. 1.2l requires no approx- 
imations or assumptions the expression for the FCNC's are exact. However, the perturbative 
approach in Section 12.1.11 provides expressions for the FCNCs that are only valid up to a 
certain order in the off-diagonal CKM matrix elements. For the perturbative approach in 
Section EHm to be valid we need the two expression for the FCNC's to be equal to at least 
quadratic order in the off-diagonal CKM matrix elements. However, as discussed above, 
comparing the results of Ref. ^H] and Ref. |2D] this is clearly not true for the (2, 1) and (1, 2) 



components of the down quark-Higgs couplings Xr^. 

In the uniform squark limit, the flavor violating coupling given in Eq. ()25|1 has the form 



^ t;,(l + e3tan/5)(l + eotan/?) ^ ' 

which does not agree with the results in Ref. [TH], where they find the corrected coupling to 
be 

^^^^^ - V, |l + 6otan/5P(l + eotan/3)2^^//^^//^'^" x.tan/:/). (28) 

To understand this difference between the results of Ref. JH] and Ref. we need to look 
at the approximations made in Ref. Diagonalizing the tree level quark mass matrices in 
Eq. ((7j) leads to uncorrected diagonal masses md and a CKM matrix Vq. However the large 
tan f] enhanced radiative corrections lead to off-diagonal terms in the mass matrix, which 
have the form 

(md + Amd)-"' = rrid, {{1 + ej tan I3)6-^' + eYyt tan pX^o') (29) 

where = Vq'^*Vq^ for J 7^ / and ej is defined in eqn. (j2Sl)- We have also neglected 
contributions to the diagonal elements of the form jV^^'^P for J 7^ 3 as they are subdominant. 
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Hence, to go to the physical quark basis we need to further diagonahze this effective mass 
matrix by unitary matrices Dl,r so that 

e-^^^(D^^(md + Amd)DL)''' = m^.d'' (30) 

where 6j = arg(l + ejtan/5). The approach taken in Ref. [IH] is to perturbatively expand 
the diagonahzation matrices Dl and Dr so that 

Dl = 1 + ADl (31) 
Dr = 1 + ADr (32) 

where the unitarity of Dl,r to hnear order in A leads to conditions (ADl,r)^ = — ADl,r, 
so that when J ^ I m Eq. ()30|) we have the condition 

e-^^-^(-(ADR)md + Amd + md(ADL))-^^ = 0, (33) 

where the ffid includes higher order terms and higher orders in A have been neglected. Using 
Eq. ()33|1 and its dagger along with the hierarchy in quark masses gives us 

(ADl)^^ = <! .j^'' . . . (34) 



and 



l+eltan/S-^O J < I 



/ A-r-v \JI I 1^dj yl+ejtan/3 L^tjuaii/jy " I'Qc;'! 

"i^dj I tan /3 l+ejtan/3 J 

Putting these matrices back into Eq. ()33|) with (J, /) = (2, 1) the dominant terms have the 
form 

I + €2 tan p 

which are comparable to the terms that were neglected in Eq. like 

e ((Amd)(AD,)) = " (1 + tan/3)(l + 63 tan^) " ^'^^ 

This is particularly true for values of €3 < and large values of tan (3. Therefore the deviation 
between Ref. [201 and Ref. [12] in the Kaon sector is due to a breakdown in the perturbative 
series leading to first and second order contributions being comparable. The expansion 
shown in Ref. [12] works for the (1, 3), (2, 3), (3, 1) and (3, 2) components as they expanded 
the mass matrices only to first order. As mentioned above, an analysis of the second order 
corrections, together with a derivation of Eqs. (fT^ -(fT ^ is presented in Appendix lA.ll 



"^rfj ( eyi/gtan/j e^ygtan/3 \ i(ej~ei)\JI J ^ J 
tan /3 y 
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2.1.4 Flavor changing in the Charged Higgs Couphng 



The process of calculating the flavor changing couplings for the charged goldstone modes is 
exactly the same as in Ref . [ISi . As the couplings of the goldstone has to match those of the 
W-bosons at tree level, so as to form its longitudinal mode, the flavor changing effects have 
to be 



72 



(38) 
(39) 



The charged Higgs has the effective lagrangian ^T] 



^eff 



V2_ 
— ur 

V 



cot iSrHu — tan /3AY„ 



VeffdLH+ + 



V2 

V 



tan /5mrf — cot /^AY^^ 

v2 



where 



(40) 



(41) 
(42) 



are the generic form of corrections to the down(up) Yukawas after neglecting the Yukawas of 
the flrst two generations. The matrices eg and e'y are closely related to cq and ey and their 
form is given in Appendix lA.il Hence, we flnd for J = 1, 2, 3 



V2 

V 



mt cot P V^lf {l - tan j3 [e^ 



l + e3tan/?^3j eyy^tan/? 
1 + eH tan /? 1 + e!] tan /? 



for J ^ 3 



Pi 

—niuj cot P V/^ 



1 - tan/5 en + ey% k ^ 

^ I + eO* tan/5 



and flnally for (J, I) = (2, 1), (1, 2), (1, 1) and (2, 2) 

V2 



(P. 



H+\JI 
RL J 



— niuj cot (3 V//f ^1 - tan f3eQ^ 



(43) 



(44) 



(45) 



which agrees with Ref. ^\ if the phases are neglected. To flnd the left-right coupling we 
neglect the ( AY^) as it is suppressed by cot P so that we have for / 7^ 3 



{PHr^' 



a/2 tan/5(l + es tan/?) /l + eo*tan/5 eyy^^tan/? 

V (1 + egtan/5)(l + e^tan/?) ^-^-^ V 1 + tan /? l + estanp 
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, (46) 



for J ^ 3 



" Vl + tan/? ^^^^ ^^^^ 



and for (J, /) = (3, 3) and J 3 / 



(pH+^33 _ \/2 mtan(3 33 

^^^^^ - Vl + 6^tan/3^^// ^^^^ 
rp//+y/ _ md,tan/3 jj 

3 Flavor changing processes in the Kaon and B^-Meson 
systems 

3.1 AF = 2 processes 

The effective Hamiltonian that contributes to AF = 2 processes in the Kaon and Bg meson 
systems have the generic form 

^e^/r=^E^^(/^)Q^(/^) (50) 

i 

where Ci{fj,) are the Wilson coefficients evaluted at the scale The AF operators for a 
meson of the form {q'^q^) are 

Qf^^ = (?M)(?i?^)- (51) 

So for the — system the quantities of interest to us are ek and the eigenstate mass 
difference AM^, which to a good approximation have the form 

AM;, = 2Re{{K'\H^ff^\K')) = -^^^^<^^'\H^fT"\K'))- (^2) 
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The SUSY contribution to the matrix element for the meson M may be written down as 



P2^^{C2^^{^SUSy) + C2^^ {fJ.su Sy)) 

+Pt''C^%liSUSY) + P^^'C^^fiSUSY)] . (53) 



For the Kaon system mx = 0.498 GeV, Fk = 0.16 GeV, the values of the NLO QCD factors 
from Ref. EOl are 



pVLL ^ Q P^^ = -18.6, P^^ = 30.6, Pf^^ = -9.3, P^^^ = -16.6 (54) 

for which the values 772 = 0.57, = 0.85 have been used. The dominant contributions 
as shown in Ref. fSl 1201 come from the double penguin diagrams which on matching give 
contributions to the Wilson coefficients 

- Gf^T^ry^Z^Mf^^^^^ ^^^^^ 

r<SRR _ °^ \^ ^ (yS \2l(vS \21 



Additional subleading contributions at large tan f3 come from the charged-Higgs boson and 
chargino box-diagram contributions to ex, and their form are given in the Appendix A. 4 of 
Ref. 

Similarly, for the Bg eigenstate mass differences AMg, using again Eq. ()50|) for AB = 2 
processes, we get, approximately, 

AM, = 2\{B,\H^/=^\Bs)\ (56) 

Therefore, using Eq. ()53|). the mass difference in the Bs — Bg meson system can be found 
using TUB, = 5.37 GeV, Fb, = 0.230 GeV and the values of NLO QCD factors from Ref. [E] 
being 

pVLL ^ g 254^ pLR ^ _o.7i, pii? ^ 0.90, Pf^^ = -0.37, P2^^^ = -0-72 (57) 

for which the values tib = 0.55, Bb^ = 1-3 have been used. Again, the dominant contributions 
come from double-penguin diagrams which have the same form as Eq. (I55j) with the indices 
(2, 1) (3, 2) and there are subdominant contributions from the box diagrams with charged 
Higgs bosons and stop-charginos. 
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3.2 AF = 1 processes contributing to Bg fi'^fi 

The effective Hamiltonian that contributes to AF = 1 processes in the Bg meson system has 
the form 

Wf/f ' = ^"•Z/'^// 1: -^'il')O.(p) (58) 



where the operators O are 



Oa 




O'a 




Os- 




O's 


= ms(bLSR)(Jl) 


Op 


= mi,{hRSL){ll->l) 


O's 





(59) 



The operators Oa and O'a can be dropped as ca and c'^ are proportional to the muon mass 
and so are small at large tan j3. Also the other primed operators are suppressed with respect 
to the unprimed ones due to the hierarchy of quark masses. So the dominant contributions 
at large tan f3 come from the penguin diagrams leading to the contributions 



Cs 



Cp 



(60) 



where O^'^ is the neutral Higgs diagonalization matrix and related to and through 
Eq. ()117|) . Hence, in the large tan/5 limit we find |^ 

Bn{Bs yuV") = 2.32 X 10~^M|^(|c5|^ + \cp\'^) (61) 



4 Numerical Results: Minimal Flavor Violation 

In this section we will study some of the phenomenological implications of the scenarios of 
minimal flavor violation. The quantities of interest in the following section are AM^, €k, 
and in particular the observables in the B sector, AMg and BlZ{Bs /i+yU~). The standard 
model theoretical prediction of AM^ has errors associated with the quantities rht,Vts, and 
Fbs \/Bbs that lead to large theoretical uncertainties (SHI EZl ■ There is good agreement 
between the central values for the SM prediction for AM^ obtained by the CKMfitter and 
UTFit groups. Their evaluation of the uncertainties is somewhat different. The UTFit group 
finds the 2cr range jHHl 

16.7 ps"^ < {AMsf^ < 26.9 ps"^ (62) 
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with central value 21.5 ps ^, which is consistent with the CKMfitter groups' 2a range 

14.9 ps-^ < (AM,)'^*^ < 31.4 ps"^ (63) 

and central value 21.7 ps"^. Additionally, the DO collaboration has reported a signal consis- 
tent with values of AM^ in the range 

21 (ps)"^ > AM, > 19 ps~^ (64) 

at the 90 % confidence level [lU]- More recently, the CDF collaboration has made a mea- 
surement of AM, j4lj, 

AM, = (17.33lgj? ± 0.07(syst))ps-\ (65) 

The experimental bound |l2l [IHl 

Bn{Bs /i+/i") < 1 X 10"^ (66) 

puts strong restrictions on possible flavor violating effects induced by the double penguin 
contributions in the large tan j3 regime. The dominant contributions for large tan (3 to AM, 
and BTZ{Bs yU^/W ) come from the same penguin diagrams. The dominant contributions to 
eg and ey come from the gluino d-squark loop and the chargino u-squark loop, respectively. 
Hence, for heavy squarks, the form of these loop corrections can be written approximately 
as 

M3||/i|Co(m?^,m|,|M3p) (67) 
AtMCo{ml,mlM^). (68) 





2a,. 


l^ol 


37r 




1 




167r2 



where Cn is the standard Passarino-Veltman function. 



4.1 Phenomenological constraints on Double Penguin Contribu- 
tions in the MFV scenario 

4.1.1 The effect of i37^(5, ) constraint on AM, 

As has been shown in Ref. [2| the chargino box diagrams can be neglected if all the squark 
masses are greater than about 0.5 TeV. We are now interested in setting an upper bound on 
the FCNC effects induced by the double penguin contributions. From the form of Eq. (jHTj) 
and Eq. ()68|) it is clear that the loop integrals are larger for smaller values of the squark 
masses. The value of eo is maximized for large values of /i and for values of M3 about twice 
the overall squark mass value. The value of ey on the other hand, is maximized for large 
values of At and values of that are of order two times the overall squark mass value. 
At the same time, large values of /x and/or may induce the presence of color breaking 
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minima |33j. Hence, values of M3 ~ 2mq ~ /i maximize ey, while pushing eo to large 
values. For these values of the parameters, the loop corrections are given by 



\i4)MAx\ ~ 2.7x10-2 (69) 
|(ey)MAx| ~ lxlO-^ (70) 

where we have constrained the trilinear mass parameter At ^ Srriq, so as not to create color 
breaking minima I45j . Let us stress that the bounds on the parameters coming from color 
breaking minima may be avoided by assuming metastability of the electroweak symmetry 
breaking vacuum. However, the somewhat extreme values of the parameters given above 
induce additional anomalies in the low energy spectrum. For instance, values of At ^ 3.2mq, 
decrease the physical Higgs mass to values lower than the current experimental bound on 
this quantity It is also important to stress that for negative values of fiM^, the 

coupling X^j^ may be enhanced by taking even larger values of \fj,\. Indeed, ey only falls off 
slowly for larger while increases linearly and therefore grows with increasing fi. 
We shall comment on the effect of taking larger values of below. 

In the region of large tan f3 the heavy CP-even and CP-odd masses are approximately 
equal and the Higgs mixing angle a ~ l/tan/3, so that the dominant contribution to 
B7l{Bs /^^/i~) is given by 

- ^^^-^ ^ (;^mS^)^ 

Similarly we find the dominant SUSY contribution to AMg comes from the C^"^ coef- 
ficient. To understand why the C^"^ term is dominant over the Cf^^ we consider the case 
when there is no CP violation in the neutral Higgs sector. In the basis (-ff°, h^, A) we have 

= (sin a, cos a, —i cos /5) and xf = (cos a, — sin a, i sin /5), where a is the Higgs mixing 
angle. Putting these values into Eq. (j33|) for the (3,2) component, we find [TUj . 

C^^ oc rhbrh, tan^ P { ""^ + + ^ ) (72) 



/ sin2(a - 


■P) 


cos2(q; - 
^ Ml. 


-P) 


V 






^sin2(Q; — 




cos2(a — 




V Ml, 




^ Mlo 





Cl^^ oc ml tan^ /5 ( ^^^^4^^ + ^2 - ^ j . (73) 

From a cursory inspection of these two equation it is not clear which term is dominant, at 
large tan /3, as Cj"-^ is suppressed by a factor of fhs/fhi, with respect to Cf^^. However, using 
the constraint equations that relate M/jO, a and j3 at tree-level in the MSSM ^1] we find 

M^^o^Mld-^) (74) 

cos2(«-/3) ^ Ml - Ml, ^ AMI .7.^ 
Ml, MliMl, - Ml,) ^ M\tan^ /3 ^ ' 
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Figure 2: Correlation between BTZ{Bs f^^f^~) and AM^. The squark masses are all 
uniform and have been set to 2 TeV. The rest of the SUSY parameters have been chosen so 
that |eo| and |ey| have their maximal values. The black lines have fixed values of Ma/ tan/?, 
but varying gluino phase. The contours represent AMg for different ranges of (M4 > 500, 
1000, 2000 GeV) for gluino mass and At phases equal to vr, and varying tan/3 values. The 
red (grey) vertical line is the experimental bound on BTZ{Bs n^fi~). The horizontal black 
line is the 2a upper bound on the double penguin contributions to AMg from the UTFit 
group while red (grey) horizontal line is the same bound from the CKMfitter group. 



where only the lowest order terms in (M|/M^) have been kept. Using these tree- level 
approximations we find that 



2 

2 



C^^ oc rhbrhs tan^ (3-^ (76) 



Cf^ocm,^tan^/3^. (77) 

Thus at large tan /? and moderate or large M^, Cl'^ clearly dominates over Cf^^? The value 
of AMs, including the corrections from new physics, may be represented as (AMs)'^^^|l + /s|, 
where /s is the total SUSY contribution. Due to C^'^ being dominant we find 

1 c 2 p2 o 

In the limit of universal squark masses, for fixed values of the supersymmetry break- 
ing mass parameters, the ratio between AMg and BTZ{Bs —>■ fi~^fi~) is proportional to 

^When the loop factors and phases are included the approximation for C^^^ still holds up to a factor of 
order 1. 
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{Ma/ Furthermore, fs is negative [T71 UHl dl . Therefore, unless > 2, the double 

penguin contributions to AM^ always interferes destructively with the SM contribution, at 
large tan/5. This result, showing the suppression of AM^ for enhanced BTZ{Bs —>■ fi~^fi~), 
has been known for some time and was first shown in Ref. pT| IT^ 

In Figs. El and 01 we show the correlation between AMg and BlZ{Bs —>■ yu^/i^) for different 
squark spectra and gaugino phases. In Fig. |21the black curves show the correlation between 
the double penguin contributions to AMg and BTZ{Bs —>■ yu^/i") for uniform squark masses 
~ 2 TeV. We have chosen the uniform squark masses to be ~ 2 TeV so as to ensure that for 
Ma < 1 TeV the effective Lagrangian in Eq.(^ and Eq.Q remains valid. Had we chosen 
squark masses of the order of 1 TeV, then the low-energy effective theory would break down 
for Ma close to 1 TeV, and a more detailed analysis of the e^'s momentum dependence would 
be required for these large values of Ma- Each of the black curves have different values of 
Ma/ tan/?. The contours represent the maximal values of lAM^j^^, for a given value of 
BR{Bs yU"*"//"), and for a given range of values of Ma- Due to the fact that for fixed Ma, 
the ratio of \AMs\^^ to BTZ{Bs goes like 1/ tan^ /?, in order to maximize | AM^I for 

any given value of BTZ{Bs —>■ /U^yU^) we need to minimize the value of tan/?. Inspection of 
the expressions given above shows that this may be achieved by choosing positive values of fi, 
arg(M3) = a.Tg{At) = n and maximal values of |eo| and |ey|. In order to define the contours 
we have taken the values of the loop corrections given in Eq. (|70|) . The horizontal black and 
red (grey) line corresponds to an upper bound on the largest possible contribution to AMg 
from new physics using the 2a values obtained by the UTFit and CKMfitter collaborations, 
Eq. (jU^ and Eq. respectively. In order to get a precise evaluation of this bound, 

a complete fit to the flavor violating processes within the MSSM should be performed, 
something that is beyond the scope of this paper. However, since in this region of parameters 
the only relevant new flavor violating contributions are from the double penguin diagrams, 
we can make an estimate of this bound in the following way: From Eq. (j62|) or Eq. (j63|) we 
have a 2-cr range that goes from values consistent with the experimentally measured value 
up to values much larger than the measured values. Therefore the negative double penguin 
contribution can be as large as the difference between the maximum allowed SM value and 
the smallest allowed experimental value. This leads to an upper bound on the magnitude 
of the double penguin contributions to AM^ of about ~ 10 ps^^ for the UTFit limits in 
Eq. (jnS) or ~ 14.5 ps~^ for the CKMfitter limits in Eq. From Fig. Hit is clear that, 

for CP-odd Higgs masses below 1 TeV, this bound does not lead to any further constraint 
beyond the one already obtained by the non-observation of the branching ratio of the decay 
Bs 

It is possible to enhance the value of AMg beyond what we have explored, by allowing 
values of > 2 rriq. If, for instance, we consider values of ^ 3mg, for the same value of 
B7l{Bs — > /i"*"/^^) we can enhance AMg by a factor ~ 1.5. This suggests that the contours in 
Figs. 121 and 01 are not strict upper bounds, and can be further enhanced, almost in a linear 
way, by pushing \n\/mg to larger values. However, due to the extreme values of the mass 
parameters selected in defining the contours, these are indicative of the upper bound on the 
the double penguin contributions to AM^ for a given value of BTZ{Bs — > fi^fi~) for natural 
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Figure 3: Same as Fig. |2l but for third generation soft supersymmetry breaking squark 
masses equal to 0.5 TeV and first and second generation squark masses equal to 5 TeV. 

values of the mass parameters. 

In Fig. 01 we depart from the limit of universal squark masses, by setting the third 
generation squark masses ~ 0.5 TeV while the first two generation squark masses are 5 
TeV, which leads to eg having its maximal value, but eg and eg being 100 times smaller. 
Hence, this splitting of the squark masses spoils the linear correlation between AMg and 
BlZ{Bs due to the different parametric dependences of and X"^^ for split 

masses. In both Figs. |2l and El the vertical red (grey) line is the experimental bound on 
Bn{Bs n+n') in Eq. 

Figs. El and El suggest that large double penguin contributions to |AMs| may not be 
obtained, for values of eg and ey close to their maximal values in Eqs. (jU^ and (f7n|). without 
violating the BTZ{Bs l^^l^ ) bound. Due to these bounds, for values of Ma < 1 TeV, the 
double penguin corrections to AM^ are restricted to be negative and relatively small, so that 
I AM, I SUSY < 4 X 10-12 QgY^ equivalently |AM,|Susy < 6 ps"!. 

The BTZ{Bs f^^f^~) bound also constrains contributions to AM^ and AM^ to values 
within experimental errors. For example, in Fig. ^ the SUSY contributions to AMk in the 
Kaon system for uniform squarks masses are below the experimental error of 6 x 10^^^ GeV 
or 0.01 ns~^, even for large values of {Ma/ tan/?)^. These results seem to be at variance 
with those obtained in Ref. pOj- This is mainly due to the fact that the authors of Ref. [20] 
represented results in regions of parameters where the value of BTZ{Bs /^^/^~) is well above 
the present limit. Observe that, to arrive at this conclusion, the new limit on BTZ{Bs 
fj'~^fi~) is essential. From Fig. |3 we can also see how the improvement in the limit on 
BTl{Bs — i> forces the double penguin contributions to |AM/^| from SUSY to be small. 

Finally, Fig. El shows similar results for ex- As happens in the case of AM^, the results for 
values of Ma < 1 TeV are far below the current experimental value of 2.282 x 10"^. 

However, within the Minimal Flavor Violation scheme, large contributions to AM, are 
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Figure 4: Correlation between BTZ{Bs — > fi^fi^) and AM^- The squark masses are all 
uniform and have been set to 2 TeV. The rest of the SUSY parameters have been chosen so 
that |eo| and |ey| have their maximal values. The black lines have fixed values of Ma/ tan/5. 
The contours are the double penguin contributions to AMk for gluino mass and At phases 
equal to tt, but varying tan/?. The left red (grey) vertical line is the present experimental 
bound on BTZ^Bg /i'^/i^) while the right blue (black) verticle line is the previous limit. 
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Figure 5: Same as Fig. HI but for ex- Only the current bound on ElZ{Bs — » /U'*"/i~) is shown, 
by the vertical red (grey) line and the horizontal blue (black) line is the experimentally 
measured value of e^. 
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Figure 6: Variation of SUSY contributions to AM,, with input parameters = 200 GeV, 
M3 = 1000 GeV, Md, = MsusY = 2000 GeV, 2Mi = M2 = M^^ = -90^ GeV^, 
At = -1000 GeV and tan [3 = 10 




Figure 7: Variation of SUSY contributions to AM^ with input parameters Mjy = 200 GeV, 
M3 = 1000 GeV, = MsusY = 2000 GeV, 2Mi = M2 = /.t, M^^ = -90^ GeV^ 

At = -1000 GeV and tan /? = 10 
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Figure 8: Variation of SUSY contributions to ex with input parameters Ma = 200 GeV, 
M3 = 1000 GeV, = MsusY = 2000 GeV, 2Mi = M2 = fi, M^^ = -90^ GeV^ 

At = -1000 GeV and tan p = 10 

possible for scenarios in which the stops and charginos are hght, so that the chargino-stop 
box diagrams become larger. Furthermore, the bound on ETZ{Bs —>■ /i"*"/^") can be satisfied 
by going to regions of large Ma or low tan (3 as chargino-stop box contributions are not very 
sensitive to tan (3. This scenario is similar to that discussed in Ref. [29^ where low values of 
tan P satisfy both the dark matter and baryogenesis constraints. In Fig. IHl we choose SUSY 
parameters 

Ma = 200 GeV, M3 = 1000 GeV, = Msusy = 2000 GeV, 
M^^ = -90^ GeV^, At = -1000 GeV tan/5 = 10, 
and 100 GeV < 2Mi, M2, fx < 500 GeV 

that agree with dark matter and baryogenesis constraints and produce a value of AMg that is 
enhanced with respect to the SM value. For this kind of particle spectrum the double penguin 
contributions to AMg are small compared to that of chargino stop diagrams. Although the 
enhancement of AMg is small, a comparison of the SM prediction, Eq. (j62|) and Eq. (j63p . 
and the experimentally measured value leads to disfavor additional positive contributions of 
AMs, larger than about 3.5 ps~^, where we have taken into account the SM allowed range 
given by the CKMfitter collaboration Eq. ()63|) . at the 2-a level. Even stronger constraints 
would be obtained if the UTfit values in Eq. ()63|) for (AMg)^^ were used. Therefore the 
smallest values of /i, smaller than 200 GeV, would be disfavored. A global fit to all flavor 
dependent observables within this scenario would be necessary in order to determine the 
precise lower bound on fi, something that is beyond the scope of this article. Also observe 
that for larger values of tan f3 there may be relevant double penguin contributions that could 
cancel the positive box-diagram contributions and therefore the bound on /i could be relaxed 
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in this case. 

Although this scenario leads to contributions to AM^ that are smaller than the present 
experimental errors on this quantity, as can be seen in Fig.[3 it leads to interesting corrections 
to exi as shown in Figure |H1 The results in Fig. IHlwere obtained for a value of the CKM 
phase (5 = 7r/3 (the best fit value within the SM). Experimentally we know that 

le^l = (2.282 ± 0.014) x 10"^ (79) 

and therefore the SUSY corrections are significant. For lower values of the CKM phase, how- 
ever, the SUSY contributions to |ei^| within this scenario can be smaller. The experimental 
value of eK is usually used to put a constraint on the p — f] plane^. The SM contributions to 
eK leads to the constraint equation jlO] 

5.3 X 10-^ = BKA^f^[{l - p)A^\%S{x*) + ^lS{xl, <) - ^{xl] (80) 

where Bk = 0.75 ±0.10, A ~ 0.85, A = 0.22, r^* = 1.32+°;^^, 77* = 0.57l[J;[}?, r^* = 0.47l|]:^^ and 
S{xt) and S{xc,Xt) are the Inami-Lim functions. Because the stops are light the dominant 
contributions to ex come from the chargino stop diagram. Under these approximations we 
find the ex constraint equation in p — f/ plane is modified to become 

5.3 X 10-^ = BkA%{1 - p)(1 + OA'XS;S{x:) + r/^S^, xl) - r/t^]. (81) 

where ( hides all the SUSY dependences. The dominant contribution to ex from SUSY 
comes from the Cvll Wilson coefficient. Thus we have approximately, 

PvLL n / 2 2 2 2 \ /oo\ 

^^8G|M^^(a:?)''^("^*V-^.V-x..-xJ ^ 

where m^-^ is the lightest stop mass, m-^^ is the lightest chargino mass and D2 is the Passarino- 
Veltmann function 

{ y\ ^ ( ^\ 

' {y — x){y — z){y — t) ^ \x) ~^ {z — x){z — y){z — t) ^\x)~^ 

log ( -) . (83) 



{t — x){t — y){t — z) \x 

Taking the lightest stop mass to be 120 GeV and approximating the the lightest chargino 
mass by \p\ we can estimate ( ~ 0.4 for values of /i ~ 100 GeV. However as IM2I ~ \fi\ there 
are also relevant contributions from the heavier chargino. Including these contributions, we 
obtain ( ~ 0.55. Including this value of ( in the theoretical prediction for ex will lead 
to a modification of the values of p and fj extracted from the fit to the fiavor observables. 
Although a global fit to these quantities within the light stop scenario is beyond the scope 
of this article, we notice that for ( < 0.55, the new constraint equation, Eq. (jHl}, is still 
consistent with the limits coming from |Kife|/|V^b|, sin(2/3)eff and AMs^d and therefore this 
scenario is not ruled out by these considerations. 

•^yO and fj are the usual corrected Wolfenstein parameters of the CKM matrix 
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4.1.2 The effect of BlZ{Bs /i^/i ) constraint on Higgs physics at the Tevatron 
and the LHC 

As shown above, in the minimal flavor violating scheme, all dominant FCNC effects at large 
tan/3 are proportional to ey, which is directly proportional to the product of the and 
At, but inversely proportional to the square of the squark masses. The FCNC effects are 
strongly enhanced for large values of tan/3 and small values of the CP-odd Higgs mass. 
The Tevatron collider is performing searches for non-standard Higgs bosons, which become 
efficient for exactly the same conditions. Therefore, in minimal flavor violating models, 
current bounds on the rate Bs impose strong constraints on the possibility of finding 

non-standard Higgs bosons at the Tevatron collider (for related studies, see Refs. |SZl^|Sni)- 
This is particularly true for large values of the At and fi parameters, for which ey is enhanced. 

Low values of the CP-odd Higgs mass are also associated with low values of the charged 
Higgs mass. These values of the charged Higgs mass induce large positive corrections to 
the branching ratio BR{b —>■ S7). Since the measured value of BR{b — > 57) agrees well 
with the SM prediction, these large charged Higgs induced corrections to the rare decay rate 
needs to be cancelled by similarly large corrections induced by supersymmetric particles. In 
minimal flavor violating schemes, these SUSY corrections are associated with stop-chargino 
loops [llj,|lHl^[SHl- For positive (negative) values of Atfi, the corrections to the amplitude 
of the decay b ^ S'-y have the same (opposite) sign to the ones associated with the charged 
Higgs corrections, and grow linearly with tan f3. Therefore, agreement of the theoretical 
predictions with the experimental values of BR{b 57) for small values of demands 
negative values of Atfi. 

Additional constraints come from the CP-even Higgs sector. For a given value of the 
overall squark masses, the mass of the lightest CP-even Higgs boson in the large tan/5 
regime depends strongly on the parameter At. In particular, this mass is maximized for a 
value of Xt = At — /i/ tan/? ~ 2.4 Msusy (where Msusy is equal to the average stop mass) 
and minimized for values of = Due to the complicated dependence of the Higgs 
boson properties on the supersymmetric mass parameters, searches for Higgs bosons at the 
Tevatron and the LHC are usually interpreted in terms of benchmark scenarios jl7j. For 
instance, the scenario with Xt/MsusY — 2.4 is named the maximal mixing scenario, since 
it is associated with the values of the stop mixing parameters that maximize the lightest 
CP-even Higgs mass. Similarly, = defines the minimal mixing scenario. While for 
the maximal mixing scenario the constraints coming from FCNC are particularly strong, no 
constraint from Bg /U^/W are expected to be obtained in the minimal mixing scenario. 

In Fig. ini we display the constraints in the M^-tan (3 plane that are induced by the 
requirement of obtaining a good agreement with the BR{b — > 57) and the non-observation 
of Bs —>■ /^'''/i^ at the Tevatron collider. The results are presented for different values of Xt 
and fi parameters and supersymmetry breaking squark masses equal to 1 TeV. The region 
of parameter space consistent with Bg fJ'^P-~ for fi = —100 GeV and fi = —200 GeV is 
below the dotted and dashed lines respectively. For each value of At, larger values of 
imply consistency with larger values of and smaller values of tan (3. On the other hand. 
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Figure 9: The dashed(dotted) line is the BlZ^Bg Ai^yU") experimental bound in the 
Ma — tan/3 plane for /x = — 200(— 100) GeV and the yellow (light grey) and blue (dark grey) 
bands are the 6 ^ 57 allowed regions for yU = —200 GeV and —100 GeV, respectively, in the 
uniform squark limit with MgusY = 1 TeV, \M^\ = 0.8 TeV, and 2Mi = M2 = 110 GeV. 
The red (grey) line is the projected CDF limit on H —>■ tt for lfb~^ luminosity. Larger 
luminosities would probe larger Ma and smaller tan (3. Also changing /i from —200 GeV to 
—100 GeV does not affect the CDF limit significantly. Figures (a),(c) and (e) have different 
values of Xt = At — /i/ tan/5 for arg(M3) = while (b), (d) and (f) have a arg(M3) = vr 
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the regions in the Ma — tan (3 plane that are consistent with the observed values of 

BR{h S7)^"P = S.SSiUis X 10"' (84) 

and the estimated theoretical uncertainty |54j 

\BR{h S7)^"^' - BR{h s-ff^\ < 1.3 x 10"^ (85) 

are given by the colored bands. For larger values of the bands move to smaller values of 
Ma or smaller values of tan (3. Actually, the approximate cancellation of the charged Higgs 
and chargino stop contributions implies a correlation between 1/M\ and At^tanjS. We have 
also plotted the projection of the CDF limit for non-standard MSSM Higgs boson inclusive 
searches in the A, H ^ tt channel for a total integrated luminosity of 1 fb~^. In order to 
obtain this limit we have used the approximate relation given in Ref. [H^ 

tan^ 3 

a{gg, bb A) x Bn{A ^ r+r") ~ a{gg, bb ^ A)sm ^ , m2^Q ' (§6) 

(1 + £3 tan + 9 

along with the Tevatron's reach for scenario of maximal mixing with /i ~ —200 GeV and a 
luminosity of 1 fb~^ shown in Ref. |i56j . 

The Tevatron collider is only sensitive to values of M^ smaller than about 300 GeV 
and values of tan (3 larger than about 40. For maximal mixing, Fig. Efa) shows that the 
constraints coming from flavor physics are sufficiently strong so as to restrict the parameter 
space consistent with the search for non-standard Higgs bosons at the Tevatron collider. On 
the other hand, for values of ~ 1 TeV, Fig. Efc) shows that one can obtain borderline 
consistency with the constraints coming from the flavor sector, but only for the smaller values 
of /i and Ma — 200 GeV. Finally, for values of At = 500 GeV or smaller. Fig. El^e) shows 
that the bounds coming from BR{b — >■ S7) are sufficiently strong as to strongly restrict the 
parameter space consistent with non-standard Higgs boson searches at the Tevatron collider. 

The situation is ameliorated for positive values of nM^, keeping negative values of ^At. 
In Figs. inib),(d) and (f) we have changed the sign of the gluino mass (the same results 
would be obtained by keeping the gluino mass fixed but changing the sign of fi and At). 
Positive values of fiM^ diminish the eg contributions and hence make the bound coming 
from B7l{Bs — > slightly less severe. The bound coming from B7l{b — > 57) is also 

improved, with the colored bands being slightly lower. Thus for < 1 TeV the region of 
Ma ~ 200 GeV, small fi and tan/5 ~ 50, that is not excluded by flavor physics, will be 
probed by the Tevatron Higgs searches in the near future. 

Finally, we consider the minimal mixing scenario, Xt ~ 0. In this case, the constraints 
coming from the non-observation of Bs become very weak, even for large values 

of \fi\. As we will explain below, this opens up an interesting possibility: The dominant 
charged Higgs contribution to the 6 57 amplitude at large tan/? is proportional to the 
charged Higgs coupling to top and bottom quarks given in Eq. (jl^ . Setting, for simplicity, 
Ab = makes the e'y ~ while 

^ ^fxM,icos'eiCo{ml,ml,Mi)+sm'9iCoiml,ml,Mi)). (87) 



25 



Therefore, in this case, the charged Higgs contribution to the B7l{b 57) becomes propor- 
tional to [sni EH 

1 - ^/iM3 tan P (cos2 ^^(^^1 , ^ , M|) + sin^ e,-Co(^L ' K ' ^3 ! 

^ ^ V—-, — ^ 

1 + es tan p 

where 6^ is the stop mixing angle. From Eq. (jHH|) we can clearly see that for large positive 
values of M^fi and tan (3, the charged Higgs amplitude can be strongly reduced. Furthermore 
when Xt ^ the chargino stop contribution to b sj is also small. Since, for these 
parameters, the BSM contributions to the BTZ{h —>■ 37) are small, the experimental limit 
in Eq. ()84j) puts only a weak constraint on the allowed value of M4. Moreover, as stressed 
above, for this parameter region Bg — >■ fi'^fi~ also provides no constraint because ~ 
implies small values of ey. 

Additionally, for the values of the parameters for which a cancellation of the charged Higgs 
contribution to BR{b —* 57) occurs, the usual bound on tan/3 that comes from requiring 
that t/b be perturbative up to the GUT scale may be relaxed: The bottom Yukawa has the 
form 

" .(l + e3tan/j) ^^'^ 

and as eatan/? is real and positive, and of order one for the cancellation to occur, the 
denominator suppresses the Yukawa for large values of tan (3. This leads to an enhancement 
of the upper bound on tan/3 coming from perturbative consistency in the bottom quark 
sector. 

In Fig. ^1 we illustrate such a scenario for different values of Because both the 
Eg fJ'^fJ'' and b s'~f constraints allow essentially any value of M4 ^ 100 GeV a large 
region of the M4 — tan (3 can be probed by the heavy MSSM Higgs searches at the Tevatron. 
Interestingly enough, the lightest Higgs boson mass is also close to the experimental bound 
rrih — 115 GeV in this region of parameters, and therefore it could be at the reach of the 
Tevatron collider searches. 

In conclusion, for minimal flavor violating schemes, the discovery of a non-standard Higgs 
signature at the Tevatron collider would point to a definite region of parameter space, with 
values of Xt of order of the squark masses or smaller. Larger values are strongly restricted 
by the present Tevatron, CLEO and B-factory experimental constraints. It is important 
to remark that, as the luminosity of the Tevatron increases, the probability of measuring 
Bg fJ'~^f^~ increases, and so does the one of measuring a non-standard Higgs boson signal. 
However, as it becomes clear from the above discussion, an improvement of the bound on 
Bg — > fJ'~^fi~ would put strong restrictions on the possibility of measuring a non-standard 
Higgs boson signature for moderate or large values of Xt. Conversely, if a Higgs boson 
signature were observed, with absence of observation of Bg yu^/i", it would imply either 
small values of X^, or a strong departure from minimal flavor violating scenarios. 

It is interesting to analyze the constraints that the non-observation of Bg — >■ fi^fi^ at 
the LHC, for a total integrated luminosity of order of 10 fb~^, would put on the MSSM 
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Figure 10: (a)-(c) Corresponds to /i = 500 1500 GeV with the blue (dark grey) band 

showing 6 — i> S7 allowed regions for these values of fi in the uniform squark mass limit with a 
common value of the squark masses MgusY = 1 TeV, M3 = 0.8 TeV, 2Mi = M2 = 110 GeV. 
The red (grey) line is the projected CDF limit on rr for lfb~^ luminosity. The dashed 
part of the projected Tevatron reach is an extrapolation of the curve, (d) Shows the effect 
of including the squark loop correction to Pj^^ vertex, proportional to eg, on 6 — > 57 rate 
for n = 1 TeV. The dashed line corresponds to the case when corrections are not included 
while the solid line corresponds to the case when they are included. 

parameter space. The projected Atlas bound on BR{Bs ^^^^^~) in this case would be of 
order 5.5 10~^ jHOl, and therefore would imply strong constraints on the My^-tan (3 parameter 
space (The final Tevatron bound, in case of non-observation of Bs — >■ /^"""/x" , assuming a total 
integrated luminosity of order 8 fb~^, will be close to 2 10~^ [6^ and therefore it will set 
similarly strong bounds on the parameter space). In order to study the possible implications 
for searches of non-standard Higgs bosons at the LHC, we have considered the projected 
reach of the CMS searches in the inclusive — $ + X, $ — r+r^ mode, at a luminosity 

of 30 fb -1 ing. 

From Fig. ^2 we can see that even for the most restrictive case of maximal mixing 
and negative values of /xMs, the bound coming from the non-observation of Bs —>■ fi~^fJ'~ 
would be consistent with the observation of a non-standard Higgs boson for small values 
of ~ 100 GeV and somewhat large values of 350 ^ M4 < 500 GeV. These bounds are 
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Figure 11: Comparison of the projected reach for non-standard Higgs bosons at the LHC 
in the inclusive pp —>■ ^ + X , ^ —>■ t^t~ mode (red (grey) hne) with the hmits that would 
be obtained in case of non-observation of the decay mode Bs — /i"*"/^" for an integrated 
luminosity of 10 fb~^ for /i = —100 GeV (dotted line) and /i = —300 GeV (dashed line). Blue 
(dark grey) and yellow (light grey) areas correspond to the bounds coming from BR{b 57) 
for /i = —100 GeV and /i = —300 GeV, respectively. The upper edge of the fi = —300 GeV 
area is denoted by the dot-dashed line. We show these results for a common value of the 
squark masses Mgusy = 1 TeV and (a) Xf = 2.4 TeV, (b) Xf = 1 TeV, and positive 
(negative) values of /iMs (fxAt), and IM3I ~ 0.8 TeV. 

strongly relaxed for smaller values of Xt. For instance, for < 1 TeV, observation of non- 
standard Higgs bosons would be still allowed for any value of Ma, provided < 300 GeV. 

5 Non-minimal Flavor Violation 
5.1 Gluino Contributions to AMg 

The results in the case of non-minimal flavor violation discussed in section 12.1.21 are quite 
similar to the case of minimal flavor violation. As in the case of MFV for large tan/3, the 
dominant contribution to AMg comes from the DP diagrams. However, in the non-minimal 
flavor violation scenario introduced here, the effects of gluino boxes can also be important 
and compete with the double penguin contributions. The appearence of the gluino-box 
contributions is a direct consequence of the quark-squark-gluino vertices not being diagonal 
in the flavor basis. In the case of uniform squarks masses these contributions disappear due 
to the CKM matrix being unitary. 

The double penguin contributions to BTZ{Bs fJ'^f^ ) in the non- minimal flavor sce- 
nario may be significantly larger than in the case of MVF. For instance, assuming that the 
third generation left-handed and right-handed down squark masses are light implies that the 
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Figure 12: Gluino box diagrams that make contributions to AM^ for the Nonminimal flavor 
violation. Diagrams (b) and (d) are possible because of the gluinos are Majorana and the 
lower diagrams have a relative sign difference with respect to the upper ones 



vertices in Eq. ()25|) are proportional to 

Xii oc V!/;V!;, + ey^ (90) 

where p = rriq^^/mq^. Therefore when the squark mass splitting is large these vertices can 
give large contributions to AM^ and BTZ{Bs — > /i"*"//"). However, the linear correlation 
between AM^ and B7l{Bs — > p'^p~) is not spoiled by the splitting of the squark masses as 
there is no flavor dependence in the factor multiplying ^djV^^jV^^j in Eq.(j2SI)- Therefore 
the BlZ{Bs fJ'^fJ''^) bound is still a severe constraint on large double penguin contributions 
to AMs like in the MFV scenario. 

An interesting case is one in which the gluino box diagrams dominate over the double 
penguin contributions to AMg for moderate values of p ~ 2 or 3. Similar to the light-stop 
scenario for MFV there are situations in which the gluino box diagram contributions are 
sizeable and the other contributions are suppressed. The double penguin contributions are 
suppressed for low values of tan/3. On the other hand, large values of /i and M2 suppress 
the stop-chargino box diagrams. Since the gluino box diagram effects are larger for small 
values of the left-handed squark and gluino masses, we shall investigate the case in which the 
third generation left-squark soft supersymmetry breaking parameters are about 100 GeV. 
To avoid the Tevatron bound on sbottoms we also assume that the lightest neutralino is 
within 20 GeV of the sbottom mass [H^- We can achieve this mass difference by choosing 
an appropriate value of Mi. For larger values of the soft SUSY breaking sbottom mass 
parameter, of about ~ 200 GeV the gluino box contribution becomes negligble. 

Light left-handed squarks tend to lead to large values of the T-parameter and hence are 
constrained by precision electroweak data. These large contributions to the T-parameter are 
induced by the large difference between the left-handed sbottom and stop masses and are 
proportional to the top quark mass. However, for some range of values of the right-handed 
stop mass parameter, these large contributions may be minimized. Indeed, for large values of 
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Figure 13: Variation of SUSY contributions to AM^ with input parameters Ma = 250 
GeV, Mst/sy = 1000 GeV, Mi = 110 GeV Ms = 1000 GeV, /i = 1100 GeV, M(^^ ^^^^^ = 
M^^£,^ = 1000 GeV, M(j^^_^^)^ = 100 GeV, A = 1110 GeV, tan/3 = 10 and all^refevant 
SUSY phases are zero, (a) Shows the variation of AM^ over small values of gluino mass, while 
(b) shows that in limit of large gluino mass we recover the SM value. 



the right handed stop mass parameter M^^ and Xt ~ M^^, the lightest stop mass becomes 
mainly left-handed and its mass is given by 



ml^Ml^+ml^-^yD\ (91) 

where is the small D-term contribution to the left-handed stop mass. Observe that for 
Xt ^ ^(jri top-quark mass contribution is strongly suppressed and hence the contribu- 
tion to the T-parameter becomes small jHSl- In our analysis we have chosen the stop mass 
parameters so that the relation Xt = M^^ is fulfilled. 

In FigHniwe see that for gluino masses below 200 GeV, the gluino-sbottom box contribu- 
tion yields a value of AM^ that is greater than the la bound coming from the SM. Similarly, 
in Fig. El we that there are large negative contributions to ex from the gluino box-diagrams 
for M3 < 200 GeV. The total value AM, drops below that of the SM, for M3 > 200 GeV, 
because of the interference between the diagrams in Fig. ^1 For the region M3 ^ 200 GeV, 
where AM, is large, the contributions to ex are also larger but negative, which seems to 
predict a total value of ex much smaller than the experimentally observed one. Therefore, 
the gluino box contributions to AM,, in this non-minimal flavor violating scenario with fla- 
vor changing effects induced by the CKM matrix elements, are generally small and are at 
most as large as those in the light stop scenario discussed above . In addition this scenario 
is in general highly contrived as the experimental constraints from light gluino and sbottom 
searches can be avoided only by going to a small corner of the MSSM parameter space. 
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Figure 14: Variation of the gluino contributions to eK with the gluino mass M3 for the same 
input parameters as in Fig. [121 

6 Conclusions 

In this article, we have studied the constraints on the parameter space of minimal flavor 
violating SUSY models coming from the latest constraints on Bg —>■ fi~^fi~, AM^, ex and 
BR{b 57). Firstly, we have shown that the analysis of the double penguin contributions 
to observables in the Kaon sector could not be done with the available formulae in the 
literature. We derived a new formula that describes well the Kaon sector contributions 
and show that the present constraints on Bg — > eliminate the possibility of inducing 

relevant double penguin corrections in this sector. Alternative contributions, coming from 
chargino and stop loop corrections can produce large contributions to ex, which, considering 
the present theoretical uncertainties, are consistent with the bounds coming from other flavor 
observables. 

We have also verified that the double penguin contributions to AMg interfere destructively 
with the SM contribution and are strongly constrained by the non-observation of Bg —>■ fi'^fi~ 
at the Tevatron collider. Analyzing the dependence of AMg on the supersymmetric loop 
corrections, we obtained upper bounds on this quantity for any given value of Bg f^^fJ'^ , 
for natural values of the supersymmetric mass parameters. We have also shown that for 
Ma < 1 TeV, under the current theoretical and experimental uncertainties, this bound 
is stronger than the bound on the new physics contributions that is obtained from the 
comparison of the SM predictions and the experimentally measured values. Finally, if the 
theoretical errors on AMg were reduced and the SM central value was to remain the same 
then negative corrections to AMg, like that of the double penguin contribution, would be 
necessary. However such double penguin corrections to AMg of about a few ps^^'s can be 
obtained only if BTZ{Bg —>■ ^ 3 x 10~^ for < 1 TeV, which is within the future 

sensitivity of the Tevatron collider. 

On the other hand, relevant, positive constributions to AMg may be obtained for light 
stops and charginos. The contributions may be as large as 25 percent of the SM values. 
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almost independently of the value of tan /3. Contrary to the double penguin contribution, 
the chargino-stop contributions are positive and they are more strongly constrained than the 
negative double penguin ones. Small values of the Higgsino mass, ^ < 200 GeV tend to be 
disfavored for mass parameters consistent with the scenario of electrowcak baryogenesis. We 
have also analysed a scenario in which there are flavor violating effects proportional to CKM 
matrix elements in the left-handed down squark-gluino vertices at tree-level. Although the 
box-diagrams may lead to significant contributions to AAf^ for sufficiently small gluino and 
down squark masses, this contributions are constrained to be small once the bounds on ex 
are taken into account. 

We have also analyzed the complementarity of these FCNC constraints with direct Teva- 
tron searches for heavy MSSM Higgs bosons. We have analyzed different scenarios and 
showed that B7l{b sj) and BlZ{Bs — > IJ''^IJ>~) puts strong constraints on the — tan/3 
plane. This study suggests that within minimal fiavor violating scenarios, the observation 
of non-standard MSSM Higgs bosons at the Tevatron collider would imply either moderate 
values of \Xt/MsusY\ ~ 1 and small values of or very small values of and large values 
of Interestingly enough, for values Xt < MgusY, the lightest CP-even Higgs boson mass 
is smaller than 120 GeV and therefore possibly at the reach of Tevatron high luminosity 
searches. 

Finally, we have analysed the implications of non-observation of Bs —>■ A*"*"/^" at the 
LHC, for a total integrated luminosity of order of 10 fb~^, on searches for non-standard 
MSSM Higgs bosons at this coUider. Even for the most restrictive case of maximal mixing 

and negative values of //M3, this situation would be consistent with the observation of a 
non-standard Higgs boson for small values of \fi\ ~ 100 GeV and somewhat large values of 
350 < Ma ^ 500 GeV. For < 1 TeV, instead, observation would be still allowed for any 
value of Ma, provided < 300 GeV. 
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A Appendix 



A.l A corrected perturbative approach for calculating FCNC 

We would like to develop a perturbative approach to calculating flavor changing vertices 
which in the limit of uniform eo should reproduce the exact result in Eq. (f?fj) . 

A. 1.1 Basic setup and notation 

As a starting point, we assume the form of the mass matrix 

(Md)-'' = m^, ((l + ejtan/3)5^^ + eyy^an/3Af). (92) 

As the off-diagonal elements are suppressed by CKM factors with respect to the diagonal 
elements we expand in terms of the CKM factors. Therefore first order terms are proportional 
to V^^*^ for J 7^ 3 and second terms are proportional to Vq^'^Vq^. Strictly speaking we should 
probably expand in the Wolfenstein parameter A and not in the CKM elements, however as 
all we want is the leading behaviour, it is sufficient to expand in terms of the CKM elements. 
So Md has both first and second order terms present and can be expanded to be 

Md = (Md)o + 5Md + md. (93) 

where S symbolizes terms linear in V^'^"^ for J ^ 3 and 6"^ symbolizes terms proportional to 
V^o*^-^Kf^ for J, / 7^ 3, so that 

(Md)^' = m,,(l + ej tan/5) (94) 

r mdjeyy! tan pVo^^* J ^ 3 = I 

(SMdY^ = I mteyy! tan pVo^^ J = 3^ I (95) 

y otherwise 

r^^M / rn,^^yy!^^^PV,''*V,'' (J, I) = (1,2), (2, 1) 

^'^^^^ = \ otherwise " ^^^> 

Now as we have second order terms explicitly in the mass matrix we need to expand the 
diagonalization matrices to second order. Additionally they have to be unitary to second 
order and the mass eigenvalues need to be real, which leads to the form 

(97) 
(98) 
(99) 
(100) 



;dl; 


= 


(1 + 5Dl + ^^Dl + 


^5Dl5Dl)^' 


;dL; 


yJI ^ 


(1 - (5Dl - (J^Dl + 


^^Dl^Dl)-'' 


Dr; 


yJI ^ 


(1 + 5Dr + ^^Dr + 


i5DR5DR)^^e^^^ 




^JI ^ 


(1 - 5Dr - ^^Dr + 


- ^5DR5DR)-^^e-*^^ 
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where ^D^ = — (5Dl,r and ^^D^ = — 5^Dl,r. Now the requirement Dl,r diagonahze 
the mass matrix Md for diagonal elements gives us the condition 

mdj ~ rud j\l + ej tan (3\ (101) 
Oj ^ arg(l + ejtan/?) (102) 

where we have only kept the leading order behaviour (i.e. 5^ terms have been neglected). 

All off-diagonal terms automatically vanish at the zeroth order and the first order con- 
tributions are the same as in Ref. 

e-e./(_(5DR)(Md)o + 5Md + (Md)o(5DL))-'^ = 0. (103) 
Which give us the results 



(*D0" = _(M^")o(.-M.)-' + (^My(M.")„ 

II. 



|(Md-^OoP-|(Md^ 



rxr, VI - (Md^")o(^M^)^^+(^Md)^^(M^- 
^'"""^ " |(Md-)oP-|(Md-)oP • 

As 5Md = for {J, I) = (1,2), (2,1) these first order corrections are zero for these ele- 
ments. To find the leading order contributions to Dl,r for these components we need to 
go to quadratic order in the expansion parameter. Therefore the condition on the leading 
contributions to Dl,r for (J, J) = (1, 2), (2, 1) are 



e-^^'^(-((52DR)(Md)o + A + (Md)o(5'DL))^' = (106) 



where 



A^^ = (5DR)-^^((5Md)=^^ + (5DL)=^^(Md)^^)-^(<5DR)-'='(5DR)=^^- 

{S^M^y^ - {6M^y'{6Bj^y' - i(5DL)''(5DL)=^^(Md)f (107) 

= i(<5DR)^=^(5DR)=^^ - {6^M^y' - {6M^y\6Bj^y' 

-^{6Bj^y'{6B^y\M^)l'. (108) 

To arrive at Eq. pOSj) we used Eq. ()103p and neglected terms of order Oim^Jmi,). Us- 
ing Eq. p06|) leads to a relation similar to the one in Eq. ()104|) and Eq. ()105p . except that 
(5Md ^ A 

a^n (M^'')o(A)^^ + (At)^^(Md^Oo 

= |(Md-)oP-|(Md-)oP 



,2r. Ui _ (M d"-^)o(At)-^^ + A^^(M ^ 
[Md^^)oP-|(Md^^^ 
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Substituting these equations into Eq. ()104|) and Eq. p05|) and neglecting all terms suppressed 
by the mass hierarchy we find 

-^ir^.vr j=3^i 

{6Bj^y'={ g|l^v^^3J* j^3 = / (111) 



otherwise 



\ tan 3 tan 3 1 / 



and 

{fhb \ 1+63 tan /3 ' l+e|tan/3 
/ ,yy2tan/3 4;/?ton/3 y i(er-fe)TA3J=^ J^n^j. (112) 
m,5 I l+ejtan/3 l+e^tan/? y / 

otherwise 

Now to calculate the leading order corrections to the (J, /) = (2, 1), (1, 2) elements we sub- 
stitute the independent and linear order terms into Eq. p09|) and Eq. ()110|) to find 



l + ea tan/3 (1 + e2tan/3)(l + egtan/?) 
ey tan^ /? 



2|l + e3tan/?|2 



(113) 



eyy/tan^ ^ e^y^^an^ ^ 



1 + 62 tan /3 1 + tan P 



|l + e3tan/5|2 (1 + tan/?)(l + tan/?) 
el^yf tan^ (3 



;i + e2tan/5)(l + e3tan;g) 



(114) 



Using Eq. ljllll) . Eq. ()112|) . Eq. pi3p and Eq. ()114p . we find the same corrections to the effective 
CKM matrix to leading order as in Refs. [191 El El E| 



ji 



e// 1+<5q tan /3 ' 



Kf/^^ ^^3 = / (115) 



V^jj otherwise 



A. 1.2 Flavor changing effective couplings of the Neutral Higgs Bosons 

Using the relations derived in the previous section, it is relatively straightforward to calculate 
the coupling of the neutral Higgs bosons to the quarks. The effective lagrangian in the initial 
basis has the form 

^eff = -{d'j)RFt^{d'j)LS' - {d'j)LK^{d'j)nS' (116) 
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where 5*° can be any of the three neutral scalars which has mixing matrix elements for the 
Higgs and xf for the Higgs. So if O^^ diagonalizes the neutral Higgs mass matrix, 
we have 

= O^^ + z sin /JO^^ 
xl = O^^ -i cos (30^^ (117) 

Now if we rotate quarks into the physical basis the Lagrangian has the form 

Ceff = -{dj)R (Dl,Ff Dl) {di)LS' - {dj)L (Dl^Ff Dk) (118) 

Therefore, assuming a mass matrix of the form given in Eq. (j92p . we obtain, 

{Yty^ = ^ {{xl + ejxl)5-^^ + eyy^X') ■ (119) 

which has a dependence up to second order on the CKM elements. Therefore, we obtain the 
following expansion in terms of CKM elements 

F^S = (Ff )o + (5Ff + (5^Ff (120) 



where 



- id J 

(Ff)o^^ = ..7-'; ,A xl + ejx'j6'^ (121) 
Vd{l + ejtanfj) 



i)£j(l+ejtan/3) 



(6Ft^y^ = ^ ^.e-S.^y?.S J = 3^1 (122) 

I i)<j(l+e3tan/3) U / 

otherwise 



^^■^Y^yi = ^ 7/iiJXn'p)^o'^*^o'* (^,/) - (1,2),(2,1) _ ^^23) 

otherwise 



Therefore the leading order contribution to the diagonal terms of ddS^ coupling is just 
Eq. p2ip . Again the zeroth term makes no contribution to the off diagonal elements of the 
ddS couplings. Hence, at linear order we have for J 7^ / 

5 (Dj,Ff Dl)^' = e~^'-' {-iSBny\Ft^Yc^ + (5Ff )JV 

{Ftyo\5B^y') (124) 

which also disappears for (J, J) = (1, 2), (2, 1). So the only contributions that are none zero 
at this order are when either J = 3 or / = 3. Using Eq. (jlllj) . Eq. ()112|) . Eq. pi5|) and 
Eq. ()122|) and neglecting terms suppressed by the mass hierarchy we find that 

{xLV = s {^iFf^^f = 



otherwise 
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where 



r 



J3 



ey (1 + el tan /?) — ey (es — ej) tan /3 



(126) 



1 + el* tan /3 

Finally to find the leading corrections to qqH coupling for (J, /) = (2, 1), (1, 2) we need 
to go to quadratic order in which case we have 



where 



r 



21 



^12 



S \21 



21 



RL) 



nisy^T^ {x^ - tan (3) 
Vdi^ + <^2 tan /?) (1 + es tan /3) 

5^ ( Dj^Ff Dl' 



^eff ^eff 



mgy'^T^^ix^ - x^ tan p) 31 32 
i;d(l + eitan/3)(l + e3tan/5) 



[(l + eotan/5)|l + e3tan/?p 



(1 + e2tan/3)|l + ei]tan/5|2 
ey?/j tan /3(1 + eg tan /3)(1 + e2 tan /5) — 
eyVt tan /?(1 + €2 tan /?)^] , 

^ — 1(1 + entan/3)|l + e3tan/3P- 

(1 + e2tan/5)|l + e[]tan/3|2 ° ^ 

eyVt tan /?(1 + €3 tan /?)(1 + €2 tan P) — eyvl tan/5(l + €2 tan /?) 



ei tan 13) + 



ei - £2 



(e^)nan2/?y2 



ey tan /? 
1 + e2 tan (3 



leyptan^ (3yl 



(127) 



(128) 



(129) 



(l + e^tan/3)(l + e^tan/3) |l + e3tan/3|2 JJ' ^^'^^^ 

In the limit that eg 's are uniform then of leading order contributions will collapse to Eq. (P7j) 
as each of F'^'^ elements go to ey. As the effective lagrangian is real, the LR couplings are 
related to RL, so that 



X 



LR 



(Xf^t 



^RL) 



(131) 



A. 2 Calculation of Loop factors 

The assumption that the squark mass matrices are block diagonal in the tree level CKM 
basis gives us 



D 



Ml 



(132) 
(133) 
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where hj = — A^j and fij = jj, — Therefore the diagonahzation matrices have the 
simple form 

_ / 5^-^ cos af^'^^ ^iJ^-^'t'f' \m.a'j^'^^ \ (^'^^\ 
^^^'^^ " 1^ sin af ^) 5" cos af '^^ ) ^'''^ 

where (pf (0^) is the phase of /t (/i) and 

cotJa? = _K'^-("1.)^ ,135) 

V2?/dj|/ij|t;« 

„t2„y = - ''"'"'^ (136) 

Following the notation of Ref. Zj^ and Z_ diagonalize the chargino mass matrix and 
Zn and diagonalize the neutralino mass matrix. Additionally, if there is a splitting in 
the mass spectrum so that the squarks of the first two generation have uniform masses (i.e. 
"^Di = = rnoi = "^Ds = rrtu^ = = uiu^ = mu^ = Msusy) we find 

4 

J^mjv, [PDC2{m\,ml^,ml^J + Q'^Co{m\,ml^,ml^J^ - (137) 
1=1 

2 ^ 

{C2{ml^, MlusY, MIusy) + MlusYC,{m\,MlusY,MlusY)) 



1=1 

2 



r2 I ,™,2 \ /,™,2 ,™,2 ,™2 \ ( 71 /f 2 ti 2 



(M4 +m^)Co(m^^,m^^3,m^^J -C2(m^„M^^5^,Mjj,5^) - (138) 



where Ci are the Passarino-Veltman functions, are the physical squark masses, Mj are 
the squark soft mass parameters and 

Pi/ = Z''^[g^Z'^-g^Zl) (139) 
Q'l = -^/ix(^-^72^I^)+^^^((M2). + my 

- ^2^^) ((M^), + {ml)j) - yliZ^fhVu. (140) 
Similarly for the antiholomorphic corrections to the up Yukawas have the form 
^o' = Y^(-^^3/^*-.C7o(|M3P,m^^,,m^,^3) + 
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J2^N, [PuC2{m%^,ml^,ml^J+Q'ilCo{ml^,mlj^,ml^j') - (141) 
1=1 

2 



1=1 

2 



-Y 



« j^^^ L yfe 
(M|,3 + ml)C^{ml^,ml^,ml^) - C2{ml^, MI^sy, MI^sy) ' 
MlusYCo{m%.Ml^SY. MI^sy)) " Z'}Z'^f,,vM<,.^l,. ^l,)] (142) 



where 



P'J = -Z'^ [giZ]^ - (143) 

g[/ = -^/ix(^+^72^^)-^^^((M^).+0 

+ 92zi^ {{Ml) J + ml) + yliZ^njv^. (144) 

The infinities present in C2 in ey's clearly cancel, however the infinities in eg's need to be 
absorbed by counter terms in the effective lagrangian. So that the C2 contributions to the 
e)'s in the above formulae are purely the finite pieces. 
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O.QOD 



J^/Tan(p)= lOGeV - 

^\/Tan(p)=20GeV 

^^/Tan(p)=30GeV - 

rv^/Tan(p)=40GeV - 

rv^/Tan((3)=50GeV - 
M^> 500 GeV □ 
rv^> 1000 GeV □ 
M^> 2000 GeV □ 



500 



5i?(B,4^|j(tp[ey|L 10 



